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We numerically study the non-equilibrium dynamics of the three dimensional Heisenberg Edwards- 
Anderson spin glass following a sudden quench to its low temperature phase. The subsequent aging 
behavior of the system is analyzed in detail, and the scaling behavior of various space-time correlation 
functions is investigated for both spin and chiral degrees of freedom. We carefully compare our 
results with those obtained from simulations of the more studied Ising version of the model, and 
with experiments on real spin glasses in which the spins have vectorial character. Finally, the 
present dynamical study offers new perspectives into the possibility of spin-chirality decoupling at 
low temperature in vectorial spin glasses. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 



I. INTRODUCTION 

Spin glass physics has been widely studied over the 
last decades because spin glasses are considered as the 
paradigm for investigating the 'glass' state**. In particu- 
lar, the Ising Edwards- Anderson spin glass model defined 
by 



Hi s 



/ , JijSiSj, 



(1) 



has been very heavily studiedi*^, since it is the simplest 
model with the necessary ingredients of randomness and 
frustration. Here, the Si = ±1 are Ising spins on a reg- 
ular lattice interacting through nearest neighbor inter- 
actions, Jij, which are random variables drawn from a 
distribution of zero mean. The poor theoretical under- 
standing of issues such as the phase diagram of the Hamil- 
tonian in Eq. jl}. the nature of its low temperature phase, 
or the extension of its mean-field solution to finite dimen- 
sions shows that the problem is indeed challenging. Also, 
due to the nature of the problem, experiments only probe 
non-equilibrium dynamics of spin glasses at low temper- 
ature because the equilibration time of a macroscopic 
sample is infinite in this region. Experiments therefore 
pertain to the field of non-equilibrium statistical mechan- 
ics^. The variety of dynamic phenomena observed in 
experiments (aging, rejuvenation, memory, etc) can be 
viewed as additional theoretical challenges'*^. 

In recent years, several theoretical approaches to 
the slow dynamics of spin glasses described the 
physics in terms of a distribution of length scales 
whose time and temperature evolution depends on 
the specific experimental protocol, leading to a good 
qualitative understanding of the dynamics of spin 
glassea 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 : 13 ! 14 ! 15 ! 16 . Early numerical studies^ 
revealed the existence of a corresponding dynamic cor- 
relation length separating small quasi-equilibrated and 



large non-equilibrated length scales. The physical rel- 
evance of these length scales was however critically dis- 
cussed only more recently, both in simulationsi£*i*i*i£*i2*i*i 
and in experimentsi2*22*2i*2£. 

However, the connection between simulations and ex- 
periments is unclear because the spins in an experimental 
spin glass have a vectorial character, so that a more nat- 
ural Hamiltonian to consider is 



H 



Si • Sj 



(2) 



where the Si are now three-component vectors of 
unit length. The Heisenberg Edwards- Anderson model 
in Eq. has been far less studied than its Ising 
counterpart. Experimentally, anisotropy induced by 
Dzyaloshinsky-Moriya interactions allows the study of 
"Ising-like" or "Heisenberg-like" samples, depending on 
its strength. Recent experiments performed on Ising and 
Heisenberg samples revealed that the distinction indeed 
matters22*22*2*V For instance, different samples behave 
quite differently even if similar temperature protocols are 
usedSSiSi. This emphasizes the need for large scale stud- 
ies of the non-equilibrium dynamics*^ of the Hamiltonian 
in Eq. (0). 

We have therefore performed detailed non-equilibrium 
simulations of the Heisenberg Edwards- Anderson spin 
glass model in three-dimensions. In this paper we discuss 
results obtained following the simplest, yet widely stud- 
ied, experimental protocol where the system is quenched 
at initial time from a high-temperature state to its spin 
glass phase. The result of temperature shift and cycling 
experiments, and the influence of finite cooling rates on 
the dynamics are the object of a future paper—. 

The paper is organized as follows. In Scc.^J we present 
the model and give technical details. The dynamics fol- 
lowing a quench is presented in Sec. lIIII Scaling behavior 
of dynamic functions is discussed in Sec. IIVI We give a 
summary of our results in the conclusion of the paper in 
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Sec. E| 



II. MODEL AND NUMERICAL DETAILS 

We numerically study the model defined by the Hamil- 
tonian in Eq. J2J) , in which the Heisenberg spins lie on the 
sites of a three-dimensional cubic lattice with N = L 3 
sites periodic boundary conditions. The random cou- 
plings are drawn from a Gaussian distribution of zero 
mean and standard deviation unity. We use a heat-bath 
algorithm — in which the updated spin has the correct 
Boltzmann distribution for the instantaneous local field. 
This method has the advantage that a change in the 
spin orientation is always made. Times will be given 
in Monte Carlo sweeps, where one Monte Carlo sweep 
represents N = L 3 spin updates. We use a rather large 
simulation box of linear size L = 60, and discuss be- 
low in more detail this choice for L. We study several 
temperatures T = 0.16, 0.15, 0.14, 0.12, 0.10, 0.08, 0.04 
and 0.02. Although all the quantities we shall study are 
self-averaging, we average over several realizations of the 
disorder, typically 15, to increase the statistics of our 
data. 

In this paper, we simulate a single type of thermal 
history, corresponding to "simple aging experiments" , as 
opposed to the increasingly complex thermal protocols 
which have been proposed in recent years^i In a sim- 
ple aging experiment, the system is initially prepared in a 
high temperature state, T T c , and suddenly quenched 
at initial time t w = in the low temperature phase, 
T <T C . The temperature is then kept constant through- 
out the experiments, while dynamical measurements are 
performed at various "waiting times" t w after the quench. 
Some "two-time" quantities correlating the system at t w 
with that a time t later are also determined. We study 
the system for a total of 10 5 sweeps (the largest value of 
t w + t) with the following 20 values of t Wl which are in a 
roughly logarithmic progression: 2, 3, 5, 9, 16, 27, 46, 80, 
139, 240, 416, 720, 1245, 2154, 3728, 6449, 11159, 19307, 
33405, 57797. 

Contrary to its Ising version, even the location and na- 
ture of the spin glass transition of the model in Eq. (0) 
have been a matter of debate, although the situation has 
clarified somewhat recently. Early simulations reported 
the existence of a zero temperature critical poin1j22i2L2& ) 
in plain contrast to experimental findings 1 ^. Kawamura 
proposed to resolve this discrepancy by introducing the 
spin-chirality decoupling scenario22i22i2LS, based on Vil- 
lain's ideas that non-collinear ground states might exist 
in systems with vector spins^ 3 -. More recently, several 
papers^^^SiSl contradicted this scenario and argued, 
for both the Gaussian and ± J versions of the model J5J, 
that spins and chiralities in fact order at the same critical 
temperature, T c > 0. Very recent simulations 37 involv- 
ing the most efficient tools used to study the Ising spin 
glass 3 ^ conclude that the present model is characterized 
by a phase transition at T c ~ 0.16, where both spins and 



chirality simultaneously freeze. This motivates the choice 
of T = 0.16 for the upper temperature in our simulations. 
The present study has three main aims. 

(1) Comparison with experiments. As explained above, 
the model in Eq. © is best suited to describe experimen- 
tal samples where the spins also have a vectorial charac- 
ter. We want therefore to compare our numerical results 
to experimental findings. Since our work is the first large 
scale simulation of the Heisenberg spin glass in the aging 
regime, we will intentionally display a wide range of nu- 
merical data covering the whole low temperature phase 
and various observables. 

(2) Comparison with the Ising spin glass. As men- 
tioned above, experiments have revealed quantitative dif- 
ferences between Heisenberg and Ising samples. There- 
fore, we also want to compare our results with the many 
aging studies of the Ising EA model, Eq. Q . 

(3) The nature of the transition in the Heisenberg 
model. Since some of the numerical works that support 
the spin-chirality decoupling scenario are performed in 
a dynamical context, we shall discuss this issue in de- 
tail, investigating the gradual freezing with time of both 
spin and chiral degrees of freedom. Following Ref.l29l we 
define chirality as 

X <t = Si+e, • (Si X Si-eJ, (3) 

where /i G {x, y, z}, and e M is a unit vector in the direc- 
tion /i. 

III. AGING DYNAMICS 

In this section, we define and study the behavior of 
various dynamical quantities that are measured during 
the aging of the system. Their scaling properties are 
discussed in detail in the following section. 

A. Energy density 

It is a central feature of glassy materials that they 
do not reach thermal equilibrium on experimental time 
scales when they are quenched to their "glassy" phase. 
The main consequence is that physical quantities keep 
evolving with time as the system tries to reach equilib- 
rium, which is known as "aging" , a term invented by the 
polymer glass community^. 

An obvious manifestation of this out-of-equilibrium dy- 
namics is therefore the time dependence of physical ob- 
servables. In our case, it is easy to follow the evolution 
of the energy density, 

e(t a ) = (±Hy (4) 

That the system ages at all temperatures studied here 
is indeed clear from Fig. ^ where at each temperature 
T G [0.02,0.16] the time dependence of e(t w ) is evident. 



3 




FIG. 1: The time dependence of the energy density for tem- 
peratures T = 0.16, 0.15, 0.14, 0.12, 0.10, 0.08, 0.04 and 0.02 
(from top to bottom) reveals the aging of the system in its 
low temperature phase, T < T c ~ 0.16. 



B. Two-time autocorrelation functions 

The evolution with time of physical observables implies 
that the dynamics of the system is not time translation- 
ally invariant. Early studies on polymeric glasses showed 
that two-time quantities reveal the aging behavior of the 
system much more strikingly 3 ^, so that two-time corre- 
lation or response functions are widely studied in aging 
glassy materials. 

The simplest two-time quantity that has been studied 
numerically in spin glasses is the autocorrelation function 
of the spins defined by 



c{t+t w ,t w ) =—j2( s i( t + t ^)- s i( t w))- 



(5) 



This function is represented as a function of the time 
difference t for various waiting times t w and various tem- 
peratures in Fig. |2] 

The first and main observation from Fig. [3 is that 
this two-time quantity is not a function of the time dif- 
ference only, as it would be in an equilibrium system, 
C(t + t w ,t w ) ^ C(t). This is just another way of saying 
that the system is aging, but much more information can 
be extracted from this correlator. 

In more detail, the shape of the curves shown in Fig. [21 
is similar to what is commonly observed in many mate- 
rials. For small time differences, t <C t w , the curves for 
various t w superpose, implying that the dynamics is time 
translationally invariant in this time regime, reminiscent 
of some sort of "local" or "quasi" equilibrium. This will 
be clarified below. Moreover, we find that for all T < 
0.10, the curves in this time regime are consistent with 
the existence of a "plateau". At T = 0.10, the plateau 
starts to be visible at time differences t ~ 10 3 — 10 4 only, 
but its existence becomes clearer at lower temperatures. 
In mathematical terms, this means that the Edwards- 
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IK!. 2: Autocorrelation function of the spins, Eq. @, for 
L — 60, as a function of the time difference t for various 
waiting times t w logarithmically spaced in the interval t w £ 
[2, 57797] (from left to right). The temperature is T = 0.14, 
0.10 and 0.04 (from top to bottom). 



Anderson parameter defined as 

qEA = lim lim lim C(t + t w ,t w ), (6) 

is finite and positive, ^ea > 0. It is important that the 
system is large enough that it can be considered effec- 
tively infinite for the values of t and t w used. If the sys- 
tem is not large enough, the overall direction of the spins 
will wander randomly during the simulation, since there 
is no energy cost to a global rotation, with the result that 
lim^oo lim^^oo C(t + t w , t w ) = at fixed L. It is clear 
from Fig.[21that ^ea is a decreasing function of tempera- 
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FIG. 3: Autocorrelation function of the spins for T — 0.05 
and L = 16 show that the use of too small a system size 
yields stationary data for t w > 10 4 , not observed in the larger 
system used in Fig. [5] 



ture, with qEA(T — ► 0) = 1. At higher temperatures, we 
do not observe a plateau, presumably because our time 
window is too small, as can be guessed from comparing 
the curves at T = 0.14 and T = 0.10 in Fig. [2 

Although a plateau is expected, because equilibrium 
measurements indicate spin glass order below T c , we 
recall that no such plateau can be unambiguously ob- 
served in the Ising spin glass in three dimensions, where 
the short time regime is well-described by a pure power 
lawlLl&. A non-zero (?ea indicates the existence of spin 
glass phase where spins are frozen in random directions, 
and the observation of a plateau in the spin autocorre- 
lation function in Fig. is the first evidence of a stan- 
dard spin glass phase in the model in Eq. J2J using non- 
equilibrium techniques. For the Ising spin glass, analo- 
gous evidence for a transition is currently missing. Ex- 
perimentally, plateaus are also hardly visible in two-time 
correlation or response functions, but this is probably due 
to the narrow experimental time window. The existence 
of a plateau can be, however, experimentally revealed 
through a scaling analysis of two-time quantities^. 

Turning to the large time regime, t 3> t w , we observe 
that curves at various waiting times do not superpose at 
all in this regime at any temperature, fully revealing the 
aging nature of the dynamics. As in many glassy materi- 
als, it is clear that the time decay of C(t+t w ,t w ) becomes 
slower when the waiting time increases. The physical in- 
terpretation is simple: since the relaxation time of the 
sample is infinite, the only relevant time scale is the age 
of the sample t w which imposes an age-dependent relax- 
ation time: the older the sample, the slower its relaxation 
becomes. We shall discuss in the next section the scaling 
behavior of these curves. 

The long time behavior found here is qualitatively very 
different from the one reported by Kawamura^ who ar- 
gues that the spin autocorrelation function becomes sta- 
tionary at large times, even at temperatures as low as 
T = 0.05. This fact was later corrected by Matsubara et 
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FIG. 4: Autocorrelation function of the chirality, Eq. Q, as 
a function of the time difference t for various waiting times 
t w logarithmically spaced in the interval t w £ [2, 57797] (from 
left to right). The temperature is T = 0.14, 0.10 and 0.04 
(from top to bottom). 



aim* who noted that a global rotation or "drift" of the 
system could affect the dynamics, and produced curves 
similar to ours by "subtracting" by hand a global rota- 
tion of the spins. This is more simply interpreted as the 
system being too small for the order of limits in Eq. © 
to apply. Our results for an L = 16 system, plotted 
in Fig. |3| are consistent with those of Kawamura^i, and 
show stationary behavior at long times. By contrast, for 
L = 60, the spin autocorrelation function shown in Fig. El 
does not reach stationarity in the same time window even 
for temperatures as high as T = 0.16 and without sub- 
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tracting a global rotation of the spins. This shows that 
one of Kawamura's numerical arguments in favor of a 
spin-chirality decoupling^ stems from data on too small 
a size. Similar finite-size effects are most probably also 
at work in Refs. l.3lU.32l 

We now present data for the autocorrelation function 
of the chirality, 

C x (t + t w ,t w ) = 1 Y^ixt (t + t w )xt{t w )). (7) 

i 

Since the system is isotropic, C x does not depend on /i 
and we have also averaged the data over the three direc- 
tions of space. Our results are shown in Fig. 0] for the 
same parameters as for the spins. The main conclusion 
from Fig. 0] is that chiralities have essentially the same 
behavior as the spins. We observe a stationary decay at 
small t, followed by a slower, waiting time dependent de- 
cay at large times. As for the spins, the appearance of 
a plateau is clear within our time window for T < 0.10, 
indicating the existence of a non-zero Edwards- Anderson 
parameter for chirality, 

<?yEA = lim lim lim C x (t + t w ,t w ), (8) 

t^OQ t w —*OG L—>00 

which also grows when T decreases. This is expected 
since we have found above that spins freeze, which implies 
that chiralities freeze as well. 



C. Time is length 

The key problem is to understand the subtle slow 
changes that the system undergoes: what does "old" or 
"young" really mean for the sample? To answer this 
question, we turn to a spatial description of the aging 
dynamics. 

The decomposition of the decay of autocorrelation 
functions into a fast stationary process and a slow non- 
stationary one directly suggests the existence of some sort 
of local equilibrium within the sample: a spin appears 
locally equilibrated (short-time dynamics) although the 
sample as a whole is still far from equilibrium and evolves 
towards equilibrium (long-time dynamics). 

It is possible to illustrate this last statement, as was 
done in the Ising casoi£. Because of the disorder, the 
spin orientations in an equilibrium configuration are ran- 
dom, so that it is impossible to detect any domain growth 
by simply looking at the spin directions. However, two 
copies of the system, a and 6, evolving independently but 
with the same realization of the disorder, will reach corre- 
lated equilibrium configurations, so that the orientation 
of the spins in one copy can be compared with those in 
the second copy. 

It is therefore useful to define the local relative orien- 
tation of the spins as 




FIG. 5: The relative orientation of the spins in two copies of 
the system, Eq. @ , is encoded on a grey scale in a 60 x 60 x 60 
simulation box at three different waiting times t w = 2, 27 and 
57797 (from top to bottom) at temperature T = 0.04. The 
growth of a local random ordering of the spins is evident. 



In Fig. Owe present snapshots where this quantity is en- 
coded on a grey scale. Comparing three successive times, 
it becomes clear that aging is nothing but the growth 
with time of a local random ordering of the spins im- 
posed by the disorder of the HamiltonianS^i. Notice 
that the "domains" observed in Fig. [3] have highly ir- 
regular boundaries, which will influence the behavior of 
the spatial correlators discussed below. 



cos9 l (t w ) = S?(t w )-S b l (t w ). 



(9) 



Next we focus on chiral degrees of freedom, and simi- 
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FIG. 6: The sign of the chiral local overlaps in two copies of 
the system, Eq. (11UI , is encoded on a black and white scale in 
a 60 x 60 x 60 simulation box at three different waiting times 
t w — 2, 27 and 57797 (from top to bottom). The temperature 
is T = 0.04, and we have chosen /j, = x. The actual config- 
urations are the same as in Fig. but the chiral ordering, 
corresponding to the spin ordering observed there, is hardly 
visible. 



larly define a chiral local overlap as 

q^(t w )=xr(tM' b (t w ). (10) 

In Fig. we present snapshots where this quantity is 
encoded on a black and white scale, i.e. we represent the 
quantity sgn(q^) for (j, = x. Different space directions 
would give similar plots. Although a chiral ordering must 
follow the spin ordering observed in Fig. [51 this is hardly 



visible by the eye and the system appears much more 
disordered in this chiral representation. We interpret this 
as stemming from the fact that spins are actually not 
"very" correlated within the dynamic correlation length 
(see the next subsection), and so the chiralities, which 
involve three spins, are even less correlated. 



D. Four-point correlation functions 

We now go beyond qualitative pictures of black and 
white domains and measure the dynamic correlation 
length associated with the mean domain size observed 
in Figs. and El 

First, we generalize the two-site, two- replica correla- 
tion function (which is therefore a "four-point" object) 
studied in the Ising cassii to the case of Heisenberg spins 
as 

C*(r,t w ) = l^{Sf(t w ) ■ S? +r {t w ) S 4 b (U) ■ S b l+r (t w )). 

i 

(ii) 

This function measures correlations of the relative orien- 
tation of two spins separated by a distance r at time t w , 
just as the structure factor does in a pure ferromagnet. 
Note that Ci(r,t w ) is invariant under global rotation of 
the spins in either copy, and so is independent of the 
wandering of the overall spin orientation which can af- 
fect the two-time autocorrelation functions discussed in 
Sec. IIII Bl However, there will be a change of behavior in 
Ci(r,t w ) when t w is sufficiently large that the dynamic 
correlation length becomes comparable to the system size 
L, since the system then equilibrates. 

We present the space dependence of C^ir, t w ) for vari- 
ous t w and three different temperatures in Fig. [7| These 
data obviously confirm the visual impression of the snap- 
shots in Fig. UJ At a given temperature, the decay 
of C&(r, t w ) with r becomes slower at larger t w , indi- 
cating the growth with time of a dynamic correlation 
length, £(T,t w ), sometimes also referred to as a "coher- 
ence length" . Physically, this means that an "older" sys- 
tem exhibits slower dynamics because of a larger dynamic 
correlation, very much as in standard coarsening phe- 
nomenal. 

A second piece information we get from Fig. [7| is that 
the growth of £(T, t w ) is strongly dependent on tempera- 
ture, since much larger length scales can be equilibrated 
at T = 0.14 than at T = 0.04. This is expected in a dis- 
ordered system where thermal activation is likely to play 
a role, and we shall quantify this statement in the next 
section. We also note that much smaller length scales are 
reached in the same time window in the Ising spin glass, 
both in three and four dimensions, where plots similar to 
Fig. [7| typically stopiii^ at r = 5-10, instead of r = 20- 
30 used here. 

A third piece of information is that the spatial decay of 
Ca(t, t w ) is clearly not exponential, since the latter would 
correspond to straight lines in the lin-log representation 
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FIG. 7: The two-site two-replica correlator for the spins de- 
fined in Eq. 11 U as a function of r for various waiting times 
logarithmically spaced in the interval t w £ [2, 57797] (from 
left to right). The temperature is T = 0.14, 0.10 and 0.04 
(from top to bottom). Note that the range of the r axis 
changes with temperature as a result of a slower growth of 
the dynamic correlation length at lower temperature. A non- 
exponential decay is also evident from these curves. 
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adopted in Fig. Moreover, a closer inspection of the 
data shows that, as for the time decay of the autocor- 
relation functions, they can be decomposed in a short- 
distance decay, r < £(T,t w ) where the curves at various 
t w merge, and a long-distance one, r > £(T,t w ), which 
becomes slower at larger t w . This confirms the intuition 
that spins are indeed in local equilibrium on short length 
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FIG. 8: The two-site two-replica correlator for the chiralitics 
11211 as a function of r for various waiting times logarithmi- 
cally spaced in the interval t w £ [2, 57797]. The temperature 
is T = 0.14, 0.10 and 0.04 (from top to bottom). For T — 0.10 
we show both C4, x ( r \\ , tw) and Ci x (r±,t w ) which exhibit simi- 
lar behavior. Note that the range of the r axis is much smaller 
than that for the spin correlator in Fig. |7] 
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scales, but that the system as a whole is not equilibrated. 
We note that even the local equilibrium part of the decay 
seems to be non-exponential, which might be connected 
to the strongly irregular nature of the domains observed 
in Fig.0 

Finally note that the correlator in Eq. is symmet- 
ric about L/2 = 30 due to periodic boundary conditions. 
A look at Fig. [7|justifies our use of a system size L = 60 
since even for the largest waiting time and the largest 
temperature studied here, we are in the regime where 
£(T,t w ) < L/2, so that our results are not affected by 
finite size effects. 

We now turn to the chiral degrees of freedom and define 
corresponding two-site, two-replica spatial correlations of 
the chirality as 

c ix (r,t w ) = ^<xr(^)xr; a r (^kr b (^)xat(^))- 

(12) 

In the following we distinguish between two correlators 
C^ x (rii , t w ) and C4 X (r±,t w ) if r is taken in a direction 
parallel or perpendicular to /i, respectively. 

Our results for the correlators (|12|l are presented in 
Fig. [S] As for the spins, the spatial decay becomes less 
fast at larger t w , indicating a gradual random ordering 
of the chiralities. Although this ordering was not vis- 
ible on the snapshots presented in Fig. appropriate 
correlators not surprisingly perform better than the eye. 
Chiral ordering is anyway expected since ordering of the 
spins implies the one of chiralities. In agreement with the 
visual observations, however, we find that spatial correla- 
tions of chiralities are much weaker than for the spins and 
correlators are numerically indistinguishable from noise 
beyond r ~ 4. As a result, we did not attempt to per- 
form a detailed scaling analysis of spatial correlations of 
the chirality. Again, we interpret this as being due to 
the fact that chiralities are less correlated because they 
involve three spins on a length scale r = 2. 

Results for the spin and chiral two-site, two-replica cor- 
relators for L = 16, the size studied by Kawamura-S, are 
shown in Fig. [5] By comparing the top of Fig. [5] with 
Fig. [7| which is for L = 60, we see that the behavior of 
the spin function is very strongly size dependent. Phys- 
ically this is because the smaller system size comes to, 
or approaches, equilibrium on the time scale of the sim- 
ulation, whereas the larger size does not. However, the 
data for the chiral function for L — 16 in the bottom 
part of Fig. is not very different from that for L = 60 
shown in Fig. [S] This observation explains why Kawa- 
mura— found that chiral autocorrelation does not stop 
aging while spin autocorrelation does. 

Although we have not been able to precisely estimate 
the dynamic correlation length associated with chiral or- 
der, the fact that we cannot measure correlations beyond 
r ~ 4, while we can estimate spin correlations up to 
r ~ 25 shows that spin correlations are much stronger 
in the whole temperature range that we have investi- 
gated, T < 0.16. This implies that we find no tempera- 
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FIG. 9: Two-site, two-replica functions for spins (top) and 
chiralities (bottom) at T = 0.05 and L — 16. While the 
bottom figure is very similar to Fig. [§] the top figure is very 
different (note in particular the different r-range used in both 
cases) . 



ture regime below T — 0.16 where chiral order manifests 
itself unaccompanied by simultaneous spin ordering as 
a spin-chirality decoupling scenario naturally predicts. 
Hence, although the present non-equilibrium approach 
says nothing about equilibrium behavior in the thermo- 
dynamic limit, an important conclusion of this whole sec- 
tion is that, when a proper system size is used, dynamical 
studies of the Heisenberg spin glass are more simply in- 
terpreted in terms of a simultaneous phase transition at 
T c ~ 0.16 for both spin and chiral degrees of freedom. 



IV. SCALING OF DYNAMIC FUNCTIONS 

The study of several space-time correlators of the pre- 
vious section leads to the conclusion that for T < 0.16, 
spins of the Heisenberg spin glass gradually freeze with 
time in random orientations dictated by the quenched 
disorder, naturally followed by chiral degrees of freedom. 
This behavior is qualitatively similar to that of the Ising 
spin glass. In this section, we study the scaling behav- 
ior of dynamic functions defined for the spin degrees of 
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FIG. 10: Rescaled four-point functions of the spins according 
to the scaling form (1131 . and using (1141 . for temperatures 
T = 0.14, 0.10 and 0.04 (from top to bottom). Parameters 
and symbols are the same as in Fig. 



freedom, to get a quantitative description of its aging dy- 
namics, and compare the results with numerical studies 
of the Ising spin glass and with experiments. 

A. Spatial correlations 

We start our scaling analysis with the study of spa- 
tial correlations of the spins. This choice is dictated by 



TABLE I: Temperature variation of the exponent a of the 
power law in Eq. (1141 of the spatial correlations, and the 
exponent \x which occurs in the scaling of the autocorrelation 
function, see Eq. 1261 . 



T 


0.16 


0.15 


0.14 


0.12 


0.10 


0.08 


0.04 


0.02 


a 


1.1 


1.05 


1.0 


0.9 


0.8 


0.8 


0.8 


0.8 


A* 


0.97 


0.98 


0.98 


1.0 


1.01 


1.03 


1.07 


1.09 



theoretical considerations, since scaling theories of aging 
dynamics show that time correlations have natural scal- 
ing forms when expressed as a function of the dynamical 
correlation length $2i±£iiii±2ii&±i, so that its knowledge is 
of primary importance in this study. 

Following the physical discussion of the previous sec- 
tion, it is natural to suggest the following decomposition 
of Ci(r, t w ) between a "locally equilibrated" and an "ag- 
ing" part, 

C 4 (r, t w ) ~ C 4 cq(r) (paging ( j^j. j ) > ( 13 ) 

with C4agi ng (a; — > 0) ~ constant, and Caging — > oo) = 
0. As in studies of the Ising spin glassi^di, we found that 
the functional forms 

C 4cq (r)^r-^ T \ (14) 

and 

C4agin g (a;) ^ exp (-x p ) , /3 > 1, (15) 

represent the data quite well, so that a plot at fixed T 
of r a Ci(r : t w ) versus the scaling variable r/l should col- 
lapse the data for all times t w . Such scaling plots are 
indeed presented in Fig. ^3 Although the data collapse 
is quite good, small deviations can be observed in these 
scaling plots. This might be due to the fact that correc- 
tions to scaling arise at small waiting times where the 
coherence length is still small, so that a scaling regime 
defined by £ ^> 1 is not entered yet. Similar scaling plots 
were obtained for the Ising spin glass 19,41 , although on a 
more restricted spatial range. 

The temperature variation of the exponent a(T) in 
Eq. (|14fl is shown in Table [I] At very low temperatures, 
T < 0.10, it seems to be roughly constant at about 0.8. 
The same trend is also found in the Ising spin glass, al- 
though there the exponent sticks to the value 0.5. At T c 
the scaling forms in Eqs. (|13fl . and (|14f) are also expected 
to hold with a related to the anomalous exponent r\ via 
the relation 

a(T c ) = d - 2 + r). (16) 

Our estimate for 77 is therefore 

77 w 0.1. (17) 

We cannot estimate error-bars on this value since it re- 
sults from a somewhat arbitrary scaling procedure. 
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As for the Ising case, we find evidence that 

lim lim Iim C 4 (r, t w ) = 0, (18) 

r— *oo t w — >oo L— >oo 

since a > 0. Here L has to be much larger than £(T, t w ), 
otherwise the system comes to equilibrium, and we get a 
non-zero value simply because the system has spin glass 
order (see again Fig.^J). 

From the point of view of the droplet picture^ (which 
has a single ground state plus those related by global 
symmetry), Eq. (|18(l is puzzling since one expects 
linv^i-voo Ci{r,t w ) = (q 2 ) cq = g| A . The data in 
Figs. [7| and [S] are clearly inconsistent with the large r 
limits being equal to (q 2 ) C q as estimated from the au- 
tocorrelations in Figs. [3 and However, Eq. I|18(l has 
been justifiedi*^ 1 - within the replica symmetry break- 
ing (many-state, "RSB") pictur e) 19 ! 41 . The argument 3 
is that an equilibrium calculation of the correlator (II 1|) 
within a replica symmetry breaking approach predicts 
power law behavior like Eq. (|14J) with a non-zero a in 
the "zero-overlap" sector—. However, this argument is 
really for equilibrium fluctuations, and it is not obvious 
how to translate this result to the non-equilibrium situ- 
ation of interest here. In particular, a restricted average 
over the "zero-overlap" sector cannot be justified by the 
sole (trivial) observation that the global overlap is zero 
in the aging regime 3 -. 

To distinguish between the single and many states pic- 
tures, one should perform local measurements and show 
that local properties are (in)consistent with the existence 
of a single equilibrium state-— The correlator in Eq. (jl 1|> 
can do this, using the (square of the) relative spin orien- 
tation as a local physical observable. Although, as dis- 
cussed above, it is not clear to us that the argument often 
given for a power law variation within RSB theory is cor- 
rect, the droplet theory makes the clear prediction- that 
Eq. (|11H tends to a constant, g| A , for t w — » oo. Hence the 
results of our simulations, which fit the power law vari- 
ation in Eq. (|14|l with a > 0, may be closer to the RSB 
"many states" picture, at least on the time and length 
scales being probed. Theoretical studies of the two-site, 
two-replica correlation function, C4, using the Migdal- 
Kadanoff approximation in the spirit of Refs— would 
be very valuable. 



B. Length scales 

From the scaling behavior of Ci{r, t w ), we obtain a 
practical definition of the dynamical correlation length 
of the spins as the quantity leading to the best collapse 
of the data in Fig. An obvious problem is that the 
results for £ are affected by the uncertainty in the value of 
a, as is also the case in the Ising spin glass. In particular, 
this makes it impossible to estimate error-bars on our 
results for I. The evolution with time of the coherence 
length at various temperatures is presented in Fig. 1111 
where a log-log representation is used. We find that I 
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FIG. 11: Growth law of the coherence length at different 
temperatures. The data are represented by points and lines. 
The straight line is a power law fit to the data at large times 
for T = 0.16, i ~ d /z , with z » 5.0. 



grows with time, in a temperature dependent manner, as 
was anticipated above. 

From non-equilibrium critical dynamics arguments, 

one expects £ ~ tU z , where z is the critical dynamic 
exponent. An activated scaling, £ ~ (Tlni^) 1 /^, is more 
naturally expected at low temperatures, where ip is the 
so-called barrier exponent— Moreover, at a given tem- 
perature T < T c , we expect a crossover from a critical 
(power-law) regime at short times to an activated regime 
at large times, which occurs at shorter times for lower 
temperatures. 

All these theoretical expectations can successfully be 
observed in our numerical results in Fig. 1111 For T = 
0.16, a power law behavior is consistent with the data at 
large enough times, with 

2«5.0. (19) 

Here again, we did not attempt to define error-bars on 
the value of the exponent. The value z « 7.0 was found 
using the same method in the Ising spin glass 19 . It is 
the smaller value of the dynamic exponent which allows 
larger length scales to be reached in the Heisenberg model 
than in the Ising model, as was noted above. 

Clearly, at low T, the long-time behavior is not consis- 
tent with a power law and the bending of the curves is 
indeed compatible with a logarithmic growth law. Note 
that even at the lowest temperature we have simulated, 
T = 0.02, the first two decades of waiting times give a 
growth law very similar to the one obtained at T — 0.16, 
indicating that the "true" activated behavior is entered 
in the last three decades of the simulations only. This is 
consistent with the corrections to scaling observed in the 
spatial correlators, recall Fig. For this reason, our 
data do not allow a precise determination of the bar- 
rier exponent ip from a logarithmic fit to the growth 
law which contains too many free parameters. Like- 
wise, the interpolation law between critical and activated 
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regimes 



10,19.20.44 



exp 



(20) 



which was used in earlier studies, does not account 
for the crossover seen in our data. In this expression, 
£ ~ \T — T c \~ v where v is the critical exponent for the 
equilibrium correlation length, and t , t and T are mi- 
croscopic length, time and energy scales respectively. 

The conclusion is that larger time scales need to be 
studied to establish a more quantitative description of 
the growth law for £, but this is difficult because it would 
require a huge amount of computer time, as one would 
have to work with even larger system sizes. 

It is nevertheless possible to compare these results with 
the Ising spin glass. In Fig. El we replot the data of 
Fig. ^2 together with published data obtained in the 
Ising Edwards- Anderson modeliS. We use the repre- 
sentation adopted in Ref. \^ to compare similar exper- 
imental data in Heisenberg and Ising samples, and plot 
£ versus S-hi(t w ) in a lin-log scale. In this representa- 
tion, data at different times and temperatures collapse if 

zT IT 

£ ~ t w c , a growth law advocated in some early stud- 
iesi&IL^. We recognize from Fig.^lthat this functional 
form indeed reasonably accounts for the Ising data but 
not at all for the Heisenberg ones. Also, Fig. ^] clearly 
confirms that larger length scales can be studied in the 
Heisenberg model than in the Ising model in three di- 
mensions. Further, we find in the Heisenberg model clear 
evidence of an activated, logarithmic growth law, shown 
by the downward curvature in the data in Figs. 1111 and 
1121 at low temperature. This effect is not observed in the 
Ising model, see Fig. 1121 

Turning to experiments, direct measurements of the 
dynamical correlation length are obviously impossible 
since the correlator in Eq. I|ll|) can not be measured ex- 
perimentally. The least indirect method we are aware of 
was first suggested in Ref. and used more recently in 
a variety of samples^ 6 -. It consists of a dynamic estimate 
of the number of spins that relax in a coherent manner 
probed through a magnetic perturbation. Results 4 ^ on 
both Heisenberg-like and Ising-like systems were gath- 
ered in the same representation as Fig. ^| (see Fig. 4 in 
Ref. l4r|) . Our results are useful to understand the trends 
found in experiments. In particular, in Ising samples, the 
dynamical correlation length is experimentally found to 
be smaller but to grow faster with increasing time than 
in Heisenberg samples^ 6 .. We will now see that this can 
be understood from the data in Fig. U21 

The Ising spin glass has a larger dynamical exponent 
z, which means that its dynamical correlation length ini- 
tially grows more slowly, since the straight-line data at 
T = T c in Figs . ITT1 and H"2l have slope l/z. At intermedi- 
ate temperatures where experiments are usually carried 
out, the system leaves the critical regime after a crossover 
time scale t c given by 




FIG. 12: Comparison of the growth of the dynamic cor- 
relation length in Heisen berg (data from Fig. 111(1 and Ising 
(published data from Ref.Ufl) models. We use the experimen- 
tal representation of Ref. 46, plotting I versus (T /T c )\u{t m ). 
Heisenberg data are represented with closed symbols for T = 
0.16, 0.15, 0.14, 0.12, 0.10, 0.08, 0.04, 0.02 (from right to left), 
Tc ~ 0.16. Ising data are represented with open symbols for 
T = 0.95, 0.855, 0.76, 0.665, 0.57, 0.475 (from right to left), 
T c ~ 0.95. 



For the Gaussian Ising and Heisenberg Edwards- 
Anderson models in three dimensions, it is found 
thalM 



10.5, ZV | Heisenberg — 5.5, 



(22) 



tc 



\T - T r . 



(21) 



in fair agreement with the experimental values reported 
in Ref. S 10.5 for the Ising system Fe .5Mn .5TiO 3 , 
and 7 and 5 respectively for the Heisenberg systems 
CdCr1.7Ino.3S4 and Ag:Mn 2.7%. Eqs. (E2> and © 
show that Heisenberg systems leave the critical region 
at much earlier times than Ising samples. For t > t c 
there is a crossover, visible in Fig. ^] and the Heisen- 
berg data in Fig. ^1 to presumably activated behavior, 
I oc (lni) 1 /"^, where t/j is the barrier exponent. It can 
be argued^ that this crossover rationalizes the apparent 
T-dependent dynamic exponent zT/T c observed in the 
Ising model. Hence, in the experimental time window, 
Heisenberg samples lie deeper in the activated regime 
where large length scales have been reached but further 
growth is extremely slow. In the range of times available 
in the simulations, we do not see the phenomenon ob- 
served experimentally 46 that £ increases more slowly for 
Heisenberg systems than for Ising systems. However, at 
longer times, such as those in the experiments, we expect 
this would occur since the Heisenberg model is clearly 
entering the activated region where £ increases only loga- 
rithmically with t. Hence our data provides a consistent 
explanation of the results of Ref. 46. 
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C. Autocorrelation functions 

Finally, we study the scaling behavior of the spin auto- 
correlation function. This function has been much stud- 
ied in Ising spin glasses, and has recently been mea- 
sured in experiments*^ which have usually focused on 
the more easily accessible thermo-remanent magnetiza- 
tion*^. Both functions are believed to exhibit similar 
scaling behavior. 

As we did for the spatial correlations, it is natural to 
decompose C(t+t w , t w ) into an equilibrium and an aging 
part, 



C(t + t w , t w ) ~ C C q(t) C; 



h(t + t w ) 
h(t w ) 



(23) 



where h(t) is an unknown function increasing with time, 
and the scaling function C ag i ng has the following lim- 
its: C a ging(a; — ► 0) = constant and C ag ing(a; — > oo) = 0. 
Two simple expectations for the aging part are the fol- 
lowing, (i) h(t) = t leading to a simple t/t w scaling 
of the aging part of the data; (ii) h(t) — l(t), the dy- 
namical correlation length, which is the basic outcome 
of scaling theories 8 , and is also obtained in coarsening 
phenomena. It is important to notice that although pop- 
ular, t/t w is not theoretically expected to hold in spin 
glasses, as already emphasized 19 . The natural choice of 
h(t) = £(t) does not lead to t/t w scaling if the dynamics 
involves barrier activation because then £(t) grows log- 
arithmically with time, though t/t w scaling is obtained 
for critical dynamics where £(t) oc t 1 ^ . Furthermore, ex- 
actly solvable models have demonstrated the possibility 
of more general scaling forms**** that we also discuss 
below. 

In analysis of experimental data, deviations from t/t w 
scaling are often described phenomenologically by the re- 
placement 



t 

f w 



f 



(24) 



as the scaling variable in C a gi n g, where \x is an exponent 
(in practice close to unity) and / a scaling function. The 
simplest choice would be f(x) — x, but it is more com- 
mon to take 



/(*) = 
so Eq. I|24l) becomes 



1 



1 — n 



(l + .x) 1 ^-! 



t 1 

— -> u(t w ,t) = - 

tw 1 — A* 



(* + tw) 



(25) 



(26) 



since this attempts to take into account the fact that the 
real age of the sample during the measurement is t + t w 
rather than t w , see e.g. Ref. IHTl This function u(t w ,t) 
corresponds to 



hit) = exp 



1 — fj, 



(27) 
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FIG. 13: Spin autocorrelation functions rescaled using 
Eqs. CT. and I^SJ. The temperature T and the ex- 

ponent fi are indicated in each figure. 



in Eq. (j2*5)l . and reduces to t w ~^(t/t w ) for t <^t w . 

Contrary to the Ising spin glass, we find that the initial 
decay of the correlation function is better described by a 
logarithmic behavior than by a power law, and we take 



Ceq(t) = a - b log(t) = g(t), 



(28) 



where a and b are two temperature dependent parame- 
ters describes the data very well. This can be recognized 
from the initial linear aspect of the data in Fig. |3 where a 
log-lin representation was used. Experimentally also, an 
initial power law decay of the data is generally adopted, 
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but the numerical value of this exponent is extremely 
small, typically < 0.10, so that differences from a loga- 
rithmic behavior are probably small. 

Although the low temperature data in Fig. J3J) show 
evidence for a plateau, i.e. a non-zero value at large t w 
and t, Eqs. (|23|) and H28JI do not lead to a true plateau. 
However, Eq. I|23|) is sufficient to fit the data and has one 
less parameter than a function with a plateau. Also, as 
discussed many times in the Ising cas o 19 i 52 , the use of 
an additive scaling in Eq. (|2*3*|) is in principle more ap- 
propriate, but it too contains one more free parameter. 
In practice, data taken over a limited time window are 
usually not sufficient to independently fix all free param- 
eters when such an additive scaling is used££. We find 
indeed that our data can be equally well fitted with an 
additive and a multiplicative form and both lead to sim- 
ilar results as far as the scaling properties of the aging 
part are concerned. Here, we shall present results with a 
multiplicative scaling only. 

With (a, b, /i) as free parameters, we have been able 
to get a satisfactory collapse of the numerical data at all 
waiting times and temperatures by plotting g(t)C(t + 
t w ,t w ) as a function of the scaling variable u(t w , t). Rep- 
resentative results are shown in Fig. 1131 

We find that the scaling function C ag ing(aO behaves at 
large x as 

Caging (x) ~x~ Al , T<T C , (29) 

C agin g(z) ~ (logzr A2 , T < T c , (30) 

with two temperature dependent exponents Ai,A2 > 0. 
Wc find for instance that Ai(T = 0.14) « 0.5, A X (T 
0.12 = 0.4, A 2 (T = 0.04) « 2, and A 2 (T = 0.02) w 1.5. 
The existence of two different scaling forms is another 
evidence that the dynamics crosses over from a critical to 
an activated regime when T is decreased from T c down to 
0. Again, this is not the case in the Ising spin glass which 
is well described by a power law of the autocorrelation 
at large times in the whole low temperature phase, T < 

rp 17 

An outcome of the scaling procedure of Fig.^|is there- 
fore an estimate of the temperature dependence of the 
exponent /x, and we report this in Tabled We observe a 
systematic increase of fi when the temperature decreases. 
For T = T c , we find fi < 1, a behavior which has been 
called "subaging". When decreasing the temperature, 
however, fj, becomes equal to 1 around T « 0.12, and 
"superaging" , corresponding to \x > 1, is found at still 
lower temperatures. 

In the Ising spin glass, the situation is again different 
since data in three dimensions indicate that the scaling 
variable t/t w can be used, corresponding to /i = 1, in the 
whole low temperature phas o 19 i 53 . In four dimensions, 
however, the same tendency as here is found with fj, = 1 
close to T c but \i > 1 at lower temperatures^. More puz- 
zling is the experimental situation, since there one has 
fx < 1 and the temperature dependence is the opposite 
with n decreasing when the temperature is lowered^. As 
discussed in Refs. ll9u5ll however, experimental quenches 



involve a finite cooling rate, so that laboratory aging ex- 
periments are in fact temperature-shift protocols. We 
will therefore come back to subaging behavior found in 
experiments in our future paper—. 

This multiplicity of experimental and numerical be- 
haviors obviously requires some discussion. As described 
above, there is no obvious theoretical reason to expect a 
perfect t/t w scaling in spin glasses aging at low tempera- 
tures, so that deviations from this simple behavior should 
not come as a surprise. Moreover, the rescaling obtained 
in Fig. 1131 is purely phenomenological, and must be in- 
terpreted as an effective description of the data rather 
than a fundamental one. This is consistent with Kur- 
chan's work showing that persistent superaging encoded 
by Eqs. lf2*3l and (|2*7|) with jj, > 1 is strictly impossible 5 ^. 

For T — T c , however, critical scaling with I ~ t x l z nat- 
urally implies the use of £(t+t w ) / £(t w ) and therefore t/t w 
as scaling variables. Hence, it is somewhat surprising 
that the numerics indicate instead /i(T = 0.16) = 0.98. 
Since fi ~ 1 is found at a slightly lower temperature, a 
possible interpretation could be that the value of T c de- 
termined in Ref.l37lis slightly larger than the real one. An 
indication that this is possible stems from a recent work 55 
where it is argued that the finite size scaling analysis of 
Ref . |3?] systematically overestimates the critical tempera- 
ture. A second plausible explanation is that our data are 
plagued by unknown corrections to scaling in our limited 
time window at T = 0.16. 

An increasing /i when temperature is decreased can be 
rationalized as follows. It is convenient to define an ef- 
fective relaxation tim o 5 i 19 by h(t w ) / b! \t w ) , see Eq. <|2lfl) . 
When h(t) ~ log(i) this becomes t w \og{t w ), which is 
greater than t w , leading to an effective superaging behav- 
ior since the effective relaxation time grows faster than 
t w . The conclusion is that a crossover towards activated 
behavior at low T naturally implies an apparent super- 
aging behavior well consistent with our numerics. This 
argument was already used in Ref. to interpret data 
in the four dimensional Ising spin glass. 

A purely activated behavior together with a simple 
scaling of the aging part of the decorrelation would imply 
that the use of log(t + t w ) / \og(t w ) can rescale the data 
at large times 8 . We test this idea in the middle frame 
in Fig. 1141 A comparison with the top frame where t/t w 
is used shows that a logarithmic rescaling of the data is 
much better at large times. To our knowledge, this is 
the first time that a pure logarithmic rescaling of two- 
time quantities has successfully been employed in a spin 
glass simulation. This result is important because one 
of the salient predictions of the scaling approach to spin 
glasses is the prediction of activated dynamics and loga- 
rithmic scaling in the aging regime of finite dimensional 
spin glasses 8 . Experimentally, a logarithmic rescaling of 
the thermo-remanent magnetization does not collapse the 
dat a 25 i 56 . It does, however, rescale data for the out-of- 
phase susceptibility at different frequencies q i 21 i 56 , but in 
a much more restricted time regime, u)t w ^> 1. 

Theoretical insights are also provided by the exact 
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FIG. 14: Spin autocorrelation functions at T = 0.04 rescaled 
using t/t w , log(t + t w )/lag(t w ) or \og(t) / log(t w ) as scaling 
variables for the aging part (from top to bottom). We only 
display the 12 largest waiting times, i.e. t w £ [139,57797]. 



solution of the non-equilibrium dynamics of mean-field 
spin glasses which has revealed the existence of dynamic 
ultrametricitj«22i£fli5I. Physically, ultrametricity implies 
the existence of a broad distribution of relaxation times 
organized in a hierarchical manner, in the very precise 
sense described in Ref. Technically, this generalizes 
Eq. (|23|l to a continuum of scaling functions C a gi ng asso- 
ciated with a continuum of functions h(t). Interestingly, 
dynamic ultrametricity also arises in another mean-field 
solvable model which has the advantage that an exact 



form for the time decay can be worked ou1i£&. There it 
was found that two-time correlators scale with the vari- 
able log(t)/ log(t w ). This is subtly different from the log- 
arithmic rescaling used above but the two scaling vari- 
ables are not compatible, and only the latter leads to 
dynamic ultrametricity. We have applied this alterna- 
tive rescaling in the lower frame in Fig. 1141 which shows 
that it is marginally superior to the logarithmic rescaling 
used above. Note however that this alternative scaling 
variable "compresses" slightly more the numerical data, 
so that we regard both scalings as being of equally good 
quality. Note too that they are both compatible, in a 
restricted time window, with an effective description in 
terms of a superaging behavior. This is again the first 
time, to our knowledge, that aging data so directly hint 
at the possible existence of dynamic ultrametricity in a 
three-dimensional spin glass. 



V. SUMMARY AND CONCLUSION 

We have performed the first large scale numerical 
simulations of the the aging dynamics of the three- 
dimensional Heisenberg spin glass. We have measured 
and discussed in detail the behavior of several space-time 
correlators related both to spin and chiral degrees of free- 
dom. We now summarize our main results. 

We find that the non-equilibrium dynamics of spins 
and chiralities are qualitatively very similar in the range 
T £ [0.02, 0.16] that we have investigated in detail. Two- 
time autocorrelation functions show the existence of non- 
zero Edwards- Anderson parameters for spins and chiral- 
ities, and the large time behavior is qualitatively similar 
and exhibit aging behavior, C(t + t w ,t w ) ^ C(t), indicat- 
ing the simultaneous freezing of both spins and chiralities 
in this temperature window. 

Appropriate spatial correlation functions show that ag- 
ing is due to the development with time of a spin dynam- 
ical correlation length, so that aging in spin glasses can 
be thought of as a sort of coarsening, as proposed long 
ago by Fisher and Huse&. All our results are consistent 
with the growth with time of a random ordering of the 
spins imposed by the quenched disorder. The chiralities 
naturally follow the spins but display weaker correlations. 
This is because spins correlations fall off with power of 
distance in the non-equilibrium regime, so that the cor- 
relations of the chiralities, each of which involves three 
spins, probably fall off with a faster power law. 

The system's behavior can be simply accounted for 
without invoking any decoupling between spins and chi- 
ralities, and we have explicitly shown that some of the 
numerical evidence apparently supporting such a decou- 
pling suffers from strong finite size effects. 

We find that the growth law of the spin dynamical cor- 
relation length changes from a critical power law regime 
for temperatures close to T c to a slower activated regime 
at lower temperatures. This crossover naturally influ- 
ences the scaling behavior of two-time quantities, and we 



15 



have found evidence that the low temperature behavior of 
two-time autocorrelations is more naturally interpreted 
in terms of logarithmic rescalings than the standard t/t w 
scaling commonly used. In particular, predictions from 
both droplet or mean-field approaches to spin glass dy- 
namics are equally able to rationalize the time scaling of 
dynamical functions. 

Although qualitatively similar at first sight, we have 
found that the Heisenberg Edwards- Anderson model in 
three dimensions behaves quantitatively differently from 
its Ising cousin in several ways: 

1 . The Heisenberg model clearly displays an Edwards- 
Anderson parameter in the dynamics, see Eq. © 
and Fig. whereas the Ising model does not. 

2. The growth of the coherence length is not described 
by a power law but exhibits a rapid crossover to- 
wards activated dynamics, see Figs. Illlandll2l 

3. Crossover to activated dynamics is also observed 
in the scaling of two-time quantities which are not 
described by t/t w scaling in the low temperature 
phase, see Fig. El 

4. Much larger length scales are involved in the dy- 
namics at large times and low temperatures, due 
to the smaller value of the product zv of critical 
exponents, see Eq. (|22|l . 

Interestingly, the features listed above that distinguish 
the Heisenberg from the Ising Edwards- Anderson model 
also make the aging dynamics of the Heisenberg model 
much closer to the experimental situation where results 
are generally more naturally interpreted in terms of ther- 
mally activated dynamics and logarithmically slow relax- 
ations^. As far as simple aging protocols are concerned, 
we have clearly established that aging behavior has to 
be interpreted in terms of the slow growth with time of 



a spin glass dynamic correlation length which then dic- 
tates scaling behavior of physical properties. Although 
direct experimental access to t{t w ) is hard4^, recent com- 
parative studies^ of Ising and Heisenberg samples are 
in full agreement with our findings that Ising samples 
display faster growth of a smaller dynamic correlation 
length on the time scales of interest, see Fig. El How- 
ever, despite these similarities, we have noted that the 
logarithmic scalings found in this work at low temper- 
atures and large times are not necessarily observed in 
experiments because of the finite cooling rates used in ex- 
perimentsi2i£i. Finally, our observation that the Heisen- 
berg model, as opposed to the Ising version, displays a 
clear crossover towards activated dynamics is in qualita- 
tive agreement with a series of recent studies^^^^ show- 
ing that memory effects and the influence of temperature 
shifts are much stronger in Heisenberg samples. 

These conclusions, as well as our preliminary results^, 
point to the correctness of the initial intuition which mo- 
tivated this work that the present model is better-suited 
to study and understand from a microscopic viewpoint 
also more complex thermal protocols leading to further 
non-equilibrium effects. 
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